Figure S2A and S2B
p1 = DimPlot2(Bf.combined.youngBrain, reduction = 'tsne', group.by = 'integrated_snn_res.1', label = T, pt.size = 0.01)
p2 = DimPlot2(Bf.combined.youngSC, reduction = 'tsne', group.by = 'integrated_snn_res.1', label = T, pt.size = 0.01)
plot_grid(p1, p2, ncol = 2, rel_widths = c(1,1))

Figure S2C and S2D
p1 = FeaturePlot_scCustom(Bf.combined.youngBrain, reduction = 'tsne', label = F, features = c("LOC118403099", "LOC118403104", "LOC118403096", "LOC118403100"), num_columns = 1)
p2 = FeaturePlot_scCustom(Bf.combined.youngSC, reduction = 'tsne', label = F, features = c("LOC118403099", "LOC118403104", "LOC118403096", "LOC118403100"), num_columns = 1)
plot_grid(p1, p2, ncol = 2, rel_widths = c(1,1))